Resonances in Models of Spin Dependent Point 

Interactions 

C. Cacciapuoti\ R. Carlone^, R. Figari^ 



^Czech Technical University, Doppler Institute. 
QQ Bfehova 7, 11519 Prague, Czech Republic. 

O 
P^ ^Doppler Institute, 

fj Nuclear Physics Institute, Czech Academy of Sciences. 

D 25068 Res near Prague, Czech Republic. 

Q 

0\ ^Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Napoli. 

,—1 Dipartimento di Scienze Fisiche, Universita di Napoli Federico II. 

-g^ Via Cintia 80126 Napoli, Italy. 

I 

1-^ E-mail: cacciapuoti@ujf.cas.cz, carlone@ujf.cas.cz, figari@na.infn.it. 

B 

cn Abstract 

> 

■^ In dimension d = 1,2,3 we define a family of two-channel Hamiltonians obtained as 

^^ point perturbations of the generator of the free decoupled dynamics. Within the family 

*K^ we choose two Hamiltonians, Hq and if^, giving rise respectively to the unperturbed 

1^ and to the perturbed evolution. The Hamiltonian Hq does not couple the channels and 

^^ has an eigenvalue embedded in the continuous spectrum. The Hamiltonian Hg, is a small 

OO perturbation, in resolvent sense, of Hq and exhibits a small coupling between the channels. 
T^ We take advantage of the complete solvability of our model to prove with simple 

, ^ arguments that the embedded eigenvalue of Hq shifts into a resonance for H^. In dimension 

S^ three we analyze details of the time behavior of the projection onto the region of the 

^ spectrum close to the resonance. 
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1 Introduction 

Aim of the paper is to analyze a completely solvable quantum system where formation and time 
decay properties of a metastable state can be explicitly investigated. The main technical tool we 
shall make use of is the theory of point interaction Hamiltonians in a slightly generalized setting. 
The Hamiltonians we define in the paper appeared frequently, in different forms, in the chemical 



and physical literature. In [12] they are extensively utilized, without a formal classification, 
to compute scattering data for low energy inelastic scattering of electrons by atoms (see also 
j28j for more recent applications of those ideas). In the theory of neutron scattering by nuclei 
zero range potentials with spin dependent strength were often exploited (see, e.g., [22] )■ In 
the framework of the modern theory of point interaction Hamiltonians (see [2] for a thorough 
introduction) particles with spin coupled with internal degrees of freedom and/or with magnetic 
fields were investigated by various authors ([8], [18] and references therein). 
Much closer to ours are the works of C King ([2S], [22]) (whose relation with what is done here 
will be clarified in the last section), and of P. Exner [TB]. At the best of our knowledge the 
latter author was the only one to characterize multichannel point interaction Hamiltonians in 
a rigorous way. In [16] he analyzes thoroughly scattering theory for a two channel point inter- 
action in three dimensions. Using different techniques, recently made available to the theory 
of self-adjoint extensions of symmetric operators, in this paper we give an unified presentation 
in any dimension d = 1, 2, 3 of the model-Hamiltonians introduced by Exner and we add few 
results on the decay properties of metastable states. 

In a recent paper [10] we characterized a family of self-adjoint Hamiltonians generating the 
dynamics of a quantum particle interacting, via zero range forces, with an array of localized 
quantum systems with a finite number of energy levels. Preliminary results suggest that in 
those systems it is possible to investigate details of the reduced dynamics of a particle moving 
through a quantum environment. In particular in [9] we used such kind of Hamiltonians to 
examine the evolution of a quantum particle inside a tracking chamber. 

On the other hand the same family of Hamiltonians may be used to describe quantum systems 
where the role of the environment is played by the particles. Following this idea we define and 
analyze a system made up of a localized quantum bit (a model-atom or a spin) in interaction 
with a non relativistic quantum particle. We first define an unperturbed system showing a 
ground state level and an upper energy level embedded in the continuous part of the spectrum. 
We then show how a small perturbation makes the ground state move slightly in energy and the 
upper level turn into a resonance. One can verify that any characterization of resonances apply 
to our model. The analysis of resonances and their perturbations for Schrodinger operators with 
short range potentials and for point interaction Hamiltonians was given in [3], where one can 
find also a list of references about the various aspects of resonance theory exhaustive up to the 
time of its publication. More recently a renewed interest on the subject, both for Schrodinger 
and wave equation, showed up (see, e.g., [B], [21], [30], [SS])- O^^ work deals with resonances 
generated by perturbation of eigenvalues embedded in the continuum. For recent results on 
this specific aspect of resonance theory see [TT] . For theory and applications about resonances 
induced by curvature in quantum waveguides see [I3], [H], [27] and references therein. 
For the sake of simplicity we consider here only the case of a two level atom and we will 
investigate only metastables states obtained as a perturbation of embedded eigenvalues not at 
the continuous threshold. Results on the latter case will be given elsewhere. The extension 
of our results to a generic discrete array of any multilevel localized system is straightforward. 
With a suitable choice of positions, energy levels and interaction parameters it is possible to 
define Hamiltonians with a very rich spectral structure. In particular the final Hamiltonian 
can exhibit any number of metastable and isolated bound states as well as any number of 
thresholds. On the other hand any matrix-valued Schrodinger operator with smooth potential 
could be approximated by point interaction Hamiltonians of the type we define in this paper 
with minimal changes with respect to what was done in the scalar case [20] . 
The generalization to the case of a large number of non relativistic bosons interacting with the 



localized q-bit is in progress, but no modifications in the main qualitative features of the results 
presented here are expected. It may have some pedagogical interest the analogy which exists 
between this latter system with the one consisting of a non relativistic atom coupled to the 
radiation field, analyzed in a series of important papers in the last decade (see PP, [S], [22] and 
references therein). The correspondence with this much more complex quantum system can be 
established interpreting the ionization of a particle in the model we analyze in this paper as the 
creation of a photon and the localized quantum system as the atom. "Spontaneous emission" 
would correspond to the transition from a metastable state to a state where the atom gets a 
lower energy and an ionized particle is produced. A local attractive interaction between the 
atom and the particles produces a "vacuum state" well separated from the continuous part of 
the spectrum preventing any "infrared problem" in the model-system considered here. The 
case of a ground state at the continuum threshold will be treated elsewhere. 
The analogy suggests that, in spite of their simplicity, the Hamiltonians we analyze show inter- 
esting spectral characteristics considered to be typical of more elaborate and realistic systems. 



2 Basic Notation and Preliminary Results 

This section will be devoted to introduce notation and to recall few basic results about Hamil- 
tonians we will make use of in the rest of the paper. 

The system under analysis is made up of one (non-relativistic) quantum particle and one two 
level localized subsystem (q-bit or spin). 

The two level system will be described as a spin 1/2 placed in a fixed position of space, i.e. 
a unitary vector in C^. Without loss of generality we assume that it is placed in the origin. 
We denote by a*-^-' the first Pauli matrix and with Xa ^ ^ the normalized eigenvectors of the 
operator a*^^^ 

o-^^^Xa = o-Xa o- = ±1; ||Xa||c2 = 1- 

A spin state can be written as a linear superposition ax++^X-) with a, 6 G C and |ap + |6p = 1. 
The Hilbert state-space of a system made up of one particle in dimension d and one spin 1/2 is 

H:=L2(M'^)®C2 rf=l,2,3. 

A capital Greek letter will denote a vector in 7i. Any ^ G 7i admits a decomposition of the 
following form 

where the sum runs over cr = ±. In this representation the scalar product in Ti reads 

a 

Let S be the linear operator in Ti whose domain and action are defined as follows 

D{S):=C^{R\^)®€'^ rf=l,2,3 



where C(j"(M'^\0) denotes the space of infinitely differentiable functions in M'^ with compact 
support not containing the origin. 5* is a densely defined symmetric operator in 7i. A trivial 
self-adjoint extension of S is the operator 

D{H):= H'^{R'^)^C'^ d= 1,2,3 

H := -A® Ic2 + Il2 ® /Sa'-^^ /3 G M+ , 

where H^{M.'^) denotes the standard Sobolev space of functions in L^(]R'^) together with their 
first and second generalized derivatives. H describes the free evolution of a quantum particle 
and a localized spin whose reduced dynamics is generated by the Hamiltonian Pa^^\ Notice 
that no interaction between particle and spin is considered at this stage. We fixed h = 2m = 1, 
where h is the Plank's constant and m denotes the particle mass. The constant f3 has the 
dimension of an energy which, in our units, corresponds to a reciprocal square length and it 
measures half the energy gap between the spin states x+ ^i^nd X-- The action of H on vectors 
in its domain can be written as 

Its resolvent R{z) := [H — z)~'^ acts on Ti as 

R{z)^ = 5^ ( - A - z + /3a)~V<. ®Xa ^eH; ze p{H) , 



where p{H) denotes the resolvent set of H. The integral kernel of ( — A — w) , G'^{x — x'), 
is given by 
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i ^ d = 1 
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W^\^\x\) d = 2 with wgC\M+; Im(v^) >0 (1) 



e 



d = 3 



47r X 

Here Hq (r^) is the zero-th Bessel function of third kind (also known as Hankel function of the 

first kind). We recall that Hq (r^) tends to zero as \r]\ -^ oo for Im i] > and that it has a 
logarithmic singularity in zero 

TT Z vr 

where 7 is the Euler's constant. The spectrum of H is only absolutely continuouqj in particular 

app{H) = ; aess{H) = aac{H) = [-/3, 00) . 

All the self-adjoint extensions of 5* can be obtained by Krein's theory. Details of the construction 
in dimension one and three can be found in ^U\. For the sake of completeness we give in the 
following a brief account of the procedure. 



^In this paper we will adopt the classification of the spectrum given in j34] . 



Let S* denote the adjoint of S, the deficiency indices of S are defined as n± = dim[ker(S'* ± i)]. 
It is easily seen that: (ra+,n_) = (4,4) for d = 1 and {n+,n^) = (2,2) for d = 2,3. Since 
the whole family of self-adjoint extensions of S is parameterized by the family of unitary maps 
between ker(5'* + i) and ker(5'* — i), in dimension one there is a 16 real parameters family of 
self-adjoint extensions of S, while in dimensions two and three the number of free parameters is 
reduced to 4. For d ^ 4 the deficiency indices are {n+,n-) = (0,0) meaning that the operator 
S is essentially self-adjoint and its closure coincides with the free Hamiltonian H. 
Let us define 

$^ := G"-/3- ® ;^^ zeC\R, 

where G'^{x), has been defined in (fTl). For d = 1,2,3, $^ and $1 are independent solutions of 
the equation 

5*$^ = 2$^ a = ±; <!?^eD{S*); zeC\R. (2) 

In dimension two and three {$^,$1} is a basis of ker(S'* — z). In dimension one there ex- 
ist two more independent solutions of equation ([2]), namely ^l^ := (G^^^"')' ® Xa, o" = ±, 
where {G^"^'^)' is the derivative of G^~^^. In such a case a basis of the defect space is 
{$;,$-_, $^^,$^_}. 

We will make use of 2 x 2 matrices; their elements will be labeled according to the following 
notation 



M :-- 



M(+,+) M(+,_) 



Given two matrices Mi and M2 we will denote by (M1IM2) the 2x4 matrix with the first and 
second column given by the columns of Mi and third and forth column given by the columns 
of M2. 
The following theorem holds true 

Theorem 1. Let A and B be 2 x 2 matrices such that AB* = BA* and {A\B) is of maximal 
rank . Define the operator 



D{H^^) := 1^ = ^^, ® x^ G 7^1 ^ = ^^ + 5^ q„<^l ; ^' G D{H)- z G p[E 

(J u 

a' a' 
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H^^ is self-adjoint and its resolvent, R^^{z) = [H^^ — z)"^ , is given by 

a,a',a" 

Where r'^^(^) is the 2x2 matrix 

T^''iz) = BT{z) + A, 
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d = 3. 



Remark 1. For d = 2,3, theorem [T] gives all possible self-adjoint extensions of S. For d = 1, 
Haniiltonians H^^ include only extensions of 5* for which the wave function part of vectors in 
D{H^^) is continuous in the origin. To cover the whole family of self-adjoint extensions of S, 
and include Hamiltonians defined on vectors with discontinuous wave function, we should have 
taken into account the whole bases of the defect spaces of S. For the analysis of the general 
case one can refer to [TO] . 

Remark 2. It is easily seen that for Im 2; > the operator R^^{z) — R{z) is a finite rank 
operator. This implies that aess{H^^) = aess{H) = [— /?, 00) (see, e.g., ||7j). Moreover notice 
that the pure point spectrum of H is empty. In turn this means that any A G app{H^^) must 
be a polar singularity of R^^^i^z) — R{z). 

Proof. The proof of theorem o] for d = 1,3 was given, in a more general setting, in [10]. We 
give here a sketch of the proof only for d = 2. Define two linear applications A : D{S*) -^ C^ 
and A: D{S*) ^C^, A defines the "charges" q^. 



(A*), := q^ 



lim ; — 7^^ipf^{x) 



^ = J2i^.^x.eD{S*). 



\x\^o\n{\x\) 
A defines f^, i.e., the value in the origin of the regular part of the wave function. 



A direct calculation shows that the triple (C™, A,A) is a boundary value space for S (see, e.g. 
[n]), i.e. for all ^1, ^2 e D{S*) 



^1, S*^2) - {S*^l, ^2) = 5^ (A^i)^(A^2). - (A^i)^(A^2)a 



moreover A and A are surjective. Then all self-adjoint extensions of S are given by the restric- 
tions of 5** on vectors "^ satisfying 

cr' cr' 

where A and B are two 2x2 matrices satisfying AB* = BA* and (v4|i?) has maximal rank (see, 
e.g.. Theorem 3.1.6 in [21j). This proves that operators H"^^ are self-adjoint. The resolvent 
formula ^ comes directly from the resolvent formula in ^ (see also Theorem 10 in [3T]). 
Formula (|3]) is a consequence of the fact that if^^\Ef = (i?^^(2;))~V + z"^. D 



3 The Unperturbed Hamiltonian 

Our analysis will proceed as follows. Within the family of point interaction Hamiltonians 
described in the previous section we pick up an "unperturbed" Hamiltonian Ho which exhibits 
one eigenvalue embedded in the continuous spectrum. 

In the next section we then show that inside the same family one can find a sub-family of 
self-adjoint modifications of Ho indexed by a perturbation parameter e. We prove that for 
e > the embedded eigenvalue moves off the real axis and we estimate the time decay rate 
of the corresponding metastable state. For a general approach to the problem of resonances 
generated by perturbation of eigenstates embedded in the continuous spectrum see [H] . 

Definition 1. Let — oo < a < oo. The operator Hq is defined as follows 



D{Ho) ■.= \^ en 



vl/ = vl/- + ^g,$^; ^^ = J2r.^X.eDiHy, z e p{Ho) ; 
Qa = -afa , d=l; aq^ = fa, d = 2,3\ 

Ho^ := H^' + zY,Q<rK ^eZ^(^o). 



Self-adjointness of Hq comes directly from theorem [ij by taking A = I2 and B = —al2 for d = 1, 

A = ah and 5 = I2 for d = 2, 3. 

Theorem 111 gives also the explicit formula for resolvent of Hq 



where 



Roiz) = R{z) + Y. {{Vo{z))-^)^J^l,, ■ )^l 
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d=l 



d = 2 



d = 3, 



where /i = e2(T+2™)/4. 

Notice that the matrix elements of Fq have the dimension of the inverse of the square root of 

an energy in d = 1, the dimension of the square root of an energy in d = 3 whereas they are 

dimensionless in d = 2. 

The spectrum of Ho is characterized as follows 



Theorem 2. For d = 1,2,3 the essential spectrum of Hq is 

aess{Ho) = [-/3, +00) . 
The point spectrum coincides with the set of the (real) roots of equation detro(-2) = 0. 

d = 1. If ^ a < 00 the point spectrum is empty. If — cxd < a < the point spectrum consists 
of two simple eigenvalues given by 

Eq- = -(3 - — ] £"0,+ = /^ ~ X ■ 

For all —00 < a < the lowest eigenvalue, Eq_, is below the threshold of the essential 
spectrum. For a < —2y/2j3 also the second eigenvalue, -Eo,+; is below the threshold of 
the essential spectrum, in this case the point spectrum is only discrete and the essential 
spectrum is only absolutely continuous. If —2y/2f3 ^ a < the upper eigenvalue is 
embedded in the continuous spectrum, —(3 ^ £'0,+ < 13- 

d = 2. For all < fi < 00 the point spectrum consists of two simple eigenvalues given by 

1 1 

Eq = -(3 ; -E'o,+ = (3 • 

/i /i 

The lowest eigenvalue, Eq_, is always below the threshold of the essential spectrum. For 
< 2/9/i < 1 also the second eigenvalue, Eq^^, is below the threshold of the essential 
spectrum, in this case the point spectrum is only discrete and the essential spectrum is 
only absolutely continuous. If 1 ^ 2/3/i < 00 the upper eigenvalue is embedded in the 
continuous spectrum, —[3 ^ -£0,+ < 13 ■ 

(i = 3. If ^ a < 00 the point spectrum is empty. If —00 < a < the point spectrum consists 
of two simple eigenvalues given by 

Eo,- = -/3- (47ra)2 ; Eo,+ = /3 - (Anay . 

The lowest eigenvalue, -Eq,-; is always below the threshold of the essential spectrum. For 
— cxo < a < — v^2^/(47r) also the second eigenvalue, -E'o,+; is below the threshold of essential 
spectrum, in this case the point spectrum is only discrete and the essential spectrum is 
only absolutely continuous. // — \/2^/(47r) ^ a < the upper eigenvalue is embedded in 
the continuous spectrum, —f3 ^ i5o,+ < P- 

Proof. The proof of theorem 2 is straightforward. As it was stressed in remark [2] the essential 
spectrum of Hq coincides with the essential spectrum of H. The pure point spectrum of Hq 
can only be a subset of the set of the real singularities of Ro{z) — R{z), which coincides with 
the real solutions of det To{z) = 0. 

In fact, being Tq{z) diagonal, it is immediate to check that the total spectrum of Hq is the 
union of two components obtained shifting, respectively of /3 or — /?, the spectrum of a point 
interaction Hamiltonian ha of "strength" a (see p]). In fact we have 

^0 = (ha + /?) ® n+ + {ha - /3) ® n„ , 

where II-i- are the projectors on x+ and X- respectively, II-i- = {x±^')c'^X±- The dynamics 
generated by Hq describes a particle evolving in a point potential in two independent channels 
characterized by a different energy of the spin. D 

8 



Let us denote by Qo : L'^{R'^) 
defined as follows: 



the quadratic form associated to the operator Hq. Qq is 



d=l. 



D(Qo) := {^ e h| ^ = ^ ^. ® xa ; ^^ e h\r)} 






c^ = 2,3. 



D(go) := {^ G 7^1 ^ = 5^ V^. ® x<x ; V^. = V^^ + q.G 



\~a(i . 



^^G/fi(M'^), g. gC, A<-/?} 






4 The Perturbed Hamiltonian 

The perturbed dynamics is generated by the Hamiltonian H^ defined as follows 

Definition 2. Let — oo < a < cxd and < e <C |a|. The Hamiltonian H^ has the following 
domain and action 

D{H,) := JM/ G 7^1 vl/ = v]/- + J2qaK ; ^1/^ = J^^^ ® X. e DiH); z G p(4) ; 

a a 

q± = -af± -^h d=l; 
aq± + eq:f = f± d = 2,3 j 



H,^ := H^<' + zS^ q^^l ^ E D{H, 



Again self-adjointness is a direct consequence of theorem [l| Notice that H^ has been obtained 
by adding —e (resp. e) to the off diagonal terms of matrix B (resp. A) used in the definition 
of Ho for d = 1 (resp. d = 2,3). 



The resolvent of H,, is 



R,{z) = R{z) + Y, {ir,iz)r')^JK>, ■)K ze c\i 



where 



a 



r.(^) 



e 



a;2 — ^2 
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O? — £2 



d = l 



2v/^T? a^-e^l 



r.(^) 



/^ln(v/;u(z-/3)) -Z7r/2 \ 

In (7/i(;z + /5)) - m/2 
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r.(^) 
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27r 



d = 2 



7 



+ a 



\ 



Wz + p 
47r 



+ a 



rf = 3 



It is easy to verify that for all z G C\]R, there exists Eq such that for all < e < Eq 

\\R,{z) - Ro{z)\\sin,n) ^ eC (5) 

where C is a positive constanlrj Formula (tsl) implies the uniform convergence of the resolvent 
of He to the resolvent of Hq for e ^ 0. Consequently (see, e.g., [31]) the unitary group e~*^=* 
converges strongly to e~^^°^ uniformly for t in any finite interval. In this sense Hamiltonian He 
is a small perturbation of ^o- 

We also point out the operators Ho and H^ have the same form domain. In particular, let us 
denote by Q^ the quadratic form associated to the operator H^, then D{Qe) = D{Qq) and 



d=l. 



d = 2,3. 



QM = Qo[^] +4^+(o)^-(o) + V^-(o)^+(o)] 



It is quite easily seen that the difference Qs — Qo is form-bounded with respect to Qo- 
In the following theorem we prove that for e small enough Hamiltonian H^ does not have 
embedded eigenvalues even when Hq does. In particular we show that the real pole of the 
resolvent of Hq between —jS and (3 moves in a pole of Rs in the unphysical region of the 
complex energy. In fact, due to the presence of the square roots y/z — [3 and yjz + /3, the 
resolvents Rq{z) and Re{z) can be extended to operator analytic functions on the complex 
plane with two branch-cuts starting respectively from z = (3 and z = —j3. The Riemann 
surfaces then exhibit four sheets corresponding to the four couples of branches chosen for the 
square roots. We will prove that the pole belongs to the sheet defined by Im {y/z — [3) > and 
Im {^/z + P) < 0, that is, the second Riemann sheet of ^/z + (3. The proof of the theorem is 
based on a direct analysis of the singularities of the resolvent. 



^In the following we shall always denote by C a generic positive constant. The value of C can change from 
line to line. 
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Theorem 3. Assume that 

-2^/2(3 < a <0 d=l 

1 < 2/?/i < oo d = 2 (6) 

- ^/2P/{An) <a<0 d = 3 

then there exists Eq > such that for all < e < Eq the essential spectrum of H^ is only 
absolutely continuous, 

aess{He) = (TaciHe) = [-/5, +CX)) . 

Moreover for all < e < Sq the point spectrum consists of a single eigenvalue, E^_ < —(3, 
and the analytic continuation of the resolvent Re{z) through the real axis from the semi-plane 
Im ^ > has a simple pole (resonance) in z = E^^'^ where Iin(£'J'^^) < 0. For some positive 
constant C one has jiJJ^^ — -Eq.+ I ^ Ce^ 

Proof. From the compactness of Re{z) — R{z) for any z G C\M it follows that the essential 
spectrum of Hf. coincides with [—13, oo). According to theorem 6.10 in [23] if A is an embedded 
eigenvalue for H^ the projection onto the subspace corresponding to the eigenvalue A is given 
by 

Pa = s - \im{-i5) {H,-X- i5) '^ . 



We shall prove that for all e > and for all A G [— /3, cxd) the projection Px equals zero which, 
in turn, implies that H^ has no embedded eigenvalues for all e > 0. We start proving that 
Pa is zero for any A e [— /3, oo) whenever the matrix re(A) has a bounded inverse. The proof 
easily follows from the fact that there are no eigenvalues embedded in the continuos spectrum 
of Schrodinger operators with point interactions in any dimension. For the sake of completeness 
we give here a sketch of the proof (with some detail for d = ?>). 
Let ^ belong toTi, 5 > Q and A G [— /?, oo). The following inequality holds 

||(4-A-^5)-^vi/||^||i?(A + z5)vl>|| + ^|((r,(A + ^5))-i)^^^,||(<|.^-^vi/)|||$A+..||_ (7) 

The Hamiltonian H has no embedded eigenvalues, which implies that 

\im5\\Ri\ + i5)m =0. 

To study the second term on the r.h.s. of equation ([T]) we consider separately the cases A ^ ib/3 
and A = ±/9. Let us assume first A ^ ±/5. For all C £ (0, oo) and for 5 small enough the 
following inequality holds true 

C 

IIG*^^' ||L2(ud) ^^ rf=l,2,3, 

where C is a positive constant depending on d and A and the functions C^ were defined in ([I|. 
The inequality can be proved by noticing that the Fourier transform of G^ equals {\k\'^ — w)~^ 
in any dimension. For all A 7^ ±/5 and for all a we have 

0< lim5^/2||$A+i5|| ^ ^ t/ = l,2,3. (8) 

11 



Moreover it is easy to show that for any / G L (M ), d = 1, 2, 3, and ( E (0, oo) 

lim5i/2|(G^-^/)^.(^.)|=0 (9) 

which imphes that for all \1' G 7i, a, and A 7^ ±/3 

lim6^/^\{<^^-'^,^)\=0. (10) 

(5— >0 

In fact take d = 3, R > and let B[R/5^) be the open ball of radius R/5^ centered in the 
origin, with < 6 < 1. We have 



1 . T 1 /■ \fix)\ 1 r p--\-mWC.-id)\x\ 

5^|(G'^-^',/)l^(m3)| ^5-^ / ^-^dx + 5-^ / ^ \f{x)\dx 



Yi^V^^iifiu,_. + r7l' / ifr.TiP, 



^(zi) S~\\f\\LHR^) + C[ I \f{x)\'dx 



from which formula (tol) follows. 

An analogous result holds true in d = 1,2. 



From equations rt8| and (10) it follows that for all A 7^ ±/5 



hm6\{iT,iX + ^5))-')^J\{^'^-^^m\^'t''\\=0 (H) 



whenever the matrix {ri;{z)) ^ has no singularities in z = A. 
Let us analyze the cases A = ±/5. The following inequality holds 



G' l/^2™d^ ^ C—, TTT (i— 1,2,3 



1 
where C is a positive constant which depends on d. We then have 



|$f+'^KC^ ci=l,2,3. (12) 



Moreover for all / G /.^(M'^), 

\{G~\ f)\L2{Rd-j ^ C'^xz^ll/IU2(Md) c? = 2,3 

while, proceeding in the same way as it was done in the proof of formula ^ one can prove that 

^m6'/'\iG~^'j)\L2^^^ = d=l. 

o—*0 

We then have 

\{^l'-'\n ^C'^ll^ll ci = 2,3 (13) 

and 

lim5'/^|($^^-'^^)| =0 d=l (14) 

5— »0 

for all "ii en. 

Since for d = 3, and e small enough, the matrix {r^{z))^^ is regular in ^ = ±/3, inequalities 



(12) and (13) imply that equation (11) holds for A = ±/3 in d = 3. 

12 



In (i = 1,2 the statement follows from formulae (12), (13) and (14), noticing that 

lim|detr,(±/3 + Z(5)r^ = d = 2 

and 

< lim J '\^^ '-^- < oo d=l. 

<5^0 (5V2 

We are finally left to prove that (T^lz))'^ is regular for any A G [— /3, oo) which in turn means 
that the equation detr£(2;) = has no solutions in [— /5, oo). 
The equation det Tir{z) = for d = 1, 2, 3 reads 



(tia^ 



- e^) + la^z - (3) {i{a^ - e^) + 2a^/z + (3) - Ae^^^z - f^^Jz + (3 = d=l 
( In {^/fi{z-P)) - iTx/2) ( In {^ ii{z + (3)) - m/2) - {2Txef = d = 2 



(-iV^^ + 47ra)(-iV^+/3 + 47ra) -(47re)2 = d = 3 (15) 

In the following the spectral properties of H^ are analyzed for e "small enough", which will 
mean e/\a\ <^ 1. We give details only for the case d = 3. 
Let us define the function 

f[z) := ( - ^^Jz- P + 47ra) ( - i^ z + 13 + Atto) . 

Equation det V^{z) = is equivalent to f{z) = (Ane)^. We shall write z = X whenever z is real. 
If a > 0, Im /(A) < for all A ^ -/3. For a < it is easily seen that for A ^ /?, Im /(A) > 0. 
It remains to consider the case a < and —f3 < X < (3. Being 

/'(Eo,+) = ^(V2/9-(47ra)2 + 47r|a|), 

/'(£'o,+)(-2 — -E'o,+) takes real positive values only if Im z < and Re z > -E'o,+- This means 
that the solution of f{z) = (47re)^, if such a solution exists, moves in the second Riemann sheet 
for z + P in the corresponding square root. An identical argument for d = 1 proves the same 
result. For d = 2 the imaginary part always moves toward negative values while the real part 
of the resonance is greater (smaller) than £'o,+ for £'o,+ < (-£^0,+ > 0). 
To prove that in the interval —\^2(3/{Att) ^ a < there is only one eigenvalue below —(3 we 



only need to analyze equation (15) for A < —f3, for such values of A the equation is real. We 



pose a = — |a|, then the solutions of equation (15) are given by 

r2 

la 



{,/P^/iA7r\a\) - 1) {y/-l3-X/iAn\a\) - l) = ^ A < -/? 



If condition d6| is satisfied, the function on the left hand side of the equation is positive if 
and only if A < —(3 — (Ana)'^. Moreover in that interval it is a strictly decreasing function of 
A approaching +00 for A tending to —00 . Then there is one and only one eigenvalue, E^^^ 
located below —(3. Up to the first order in e"^ the explicit formula for E^^_ reads 

E,_=Eo-- , ^^^""^^ e^ + 0{e') . 

v/2/3/(47ra)2 + 1-1 

13 



We are left to prove the last statement of the theorem, i.e., that for — v2^/(47r) ^ a < 0, 
(Ts{z)) has an analytic continuation onto the second sheet of a/^H^ across the real axis 
where it has a pole. 



Let us rewrite equation (15) as 



4:71 



\a 



Alii 



\a 



e\ 



With the following position 



i 



A-K 



\a\ 






\ot\ 



the equation reads 



e 



e 



We look for fixed points of the following recurrence procedure 



Vo 



^2/? - {Ana) 
Att 



+ « 



^k 



Vk 



T]k = \a\+i 



v/2/5-[47r(a-i + |«|) 



An 



inside the ball jr^ — r^ol < Ce"^. In the definition of r^^ the square root is the analytic continuation 
of the square root of a positive number: it is defined with positive real part, while the imaginary 
part can be positive or negative according to the fact that the argument in the square root is 
in the first or second Riemann sheet respectively. It follows that |?7fc| > |a| and |^fc| < &; and 
being 



^k+l — ^k 



we also have 



Vk Vk+l 
-2 



iVk+i-Vk) 



e 



\^k+l -^k\ ^ -^\Vk+l ~Vk\ 



We use the following estimate that is a direct consequence of the definition of rjk 

v/2/5-[47r(a + |a|F ^2/3 - [471(^-1 + |«IP 



Wk+i -Vk\ 



An 



fore some positive constant C. From (16) and (17) 



An 



e^C. 



^c\ik-i 



k~l\ 



(16) 



(17) 



l^fc+l - 61 ^ -^\Tlk+l - Tlk\ ^ ^\ik - i 



a^ 



a" 



k-l\ ■ 



Summing the series we get 



leoo-^oK 



16-^0 



with 



ei-eo 



i\a\ 



^2f3 - {Anay (\a\ + ^^2f3 - {Ana) 



rS^ + 0(5') 
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Since 

-2 



e- (47r)2£2 / ^^2/3 - (47ra)^ 



^0 = — = i^is— fI 



% 2/? \' ' 47r 

the position of the resonance is 

£;r=/3- [471(^00 + |«|)]' 

=i^o.+ + ^ ^^ - ^ ^ ^^ V2/?/(4vra)2 - 1 + 0(e^) . 

Then \E1^^ — -Eo,+ | < Ce^ and the negative imaginary part of E"^^^^ means that the pole is in 
the second Riemann sheet. 

The proof of the theorem for the cases c? = 1, 2 does not differ substantially from the one given 
in the three dimensional case. To the lowest order the eigenvalue below — /? is given by 

E, _ = Eo - - , , ^ r + 0{e'') d = 1 

' 2(V8/5/a2 + l-l) 

°' In (v/2/3/i + 1) ^ ^ 

while the resonance is 



e^a^ .£" yj%i5ja^ 



16/5 ^ 2 ' 8/?/a2 

a2 + (vr/2pr + ^2, 
with a = In (^2/?/! - l) . D 






Notice that for all d = 1, 2, 3 the shift of the lowest eigenvalue i?o,- is real and negative. The 
shift of the real part of the resonance for rf = 1, 3 is always positive, while for rf = 2 it is positive 
for 1 < 2/5/i < 2 and negative for 2/5/i > 2. 

In the following we will investigate the characteristic time decay rate associated with the pres- 
ence of sharp resonances. We will refrain to give any specific definition of a metastable state and 
we will hmit ourselves to consider the decay rate of the survival probability of a spin-up state for 
an initial state whose support in energy is contained in a small spectral neighborhood around 
the position of the eigenvalue of the free Hamiltonian embedded in the continuum. Generic ini- 
tial states inside and outside the domain of H^ would have of course different decaying behaviors 
around t = 0. 

As it is well known (see, e.g., [15], [17], [19], [32] and [33]) the small time decay rate has a 
great relevance to examine the so called Zeno and anti-Zeno behavior of the initial state. In 
the following we will not give details regarding any specific initial state. We want only to stress 
that the knowledge of the generalized eigenfunctions makes available a specific formula for any 
initial state. 

Let us recall that for e small enough the essential spectrum of Hs is only absolutely continuous 
and coincides with [— /?, oo). Then for d = 1, 2, 3 and A ^ — /3 the spectral projection of H^ on 
the absolutely continuous part of the spectrum can be defined via the Stone's formula as 

Pe{\) := s - lim -^ [^^(A - i6) - Re{\ + i6)] A G [-/3, oo) . 



<5-+o+ 27ri 
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Given a state \1/ one has that 



vl/(t) 



e"*^*P,(A)^rfA. 



An exphcit formula for -Pe(A) in terms of the generahzed eigenfunctions can be derived in all 
dimensions. Here we give details only for d = 3. From a straightforward calculation it follows 
that 

where Q is the solid angle and 

I-^IA; u) = '^'^{X; u;) + Y, ^ (A; u;, 0) (r,,+ (A));/^ G^^-'" ® x.' AG [a/3, oo) . 



In the previous formula $'^(A; uj) are the generalized eigenfunctions of H, 

$'^(A; u) := r (A; ^) ® X. ; ^ (A; ^, a;) = i^ZL^e^-VA=w?= 

47r2 

A G [a (3, oo) a; G f] . 



while G;-"" and r,,+ (A) are defined by 



"° - lim G^— "+^^ 



r,+(A)= lim TJX + i5). 



The explicit expression of the generalized eigenfunctions of -ff^ reads 



l>+(A; 



U)] 



(A-/3) 

, 3 
4:712 



^^U>V>^- ^ ^^ + 



47rj 



« 



3iv^A=^M 



AW3_|^|W4±2_|^|\_^2 47r|- 



47rj 



Am 



^iV>^\-\ 



Airi 



— \a 






X- 



X^ P, Lu eVi 



X++ 



^;iX;co) 



(A + /3)i 

A 3 

Air 2 



ILD^X+fS 



'®X~ + 



47rj 



— a 



J./X+P\-\ 



/x^ 



Ani 



— a 



qM^\a\)-e2 47r|- 

47ri 



oiv^A^^M 



47rj 



— - — \a\ ] { A — \a 



Ani 



47r| 



X+ 



X^ -(3, uj en 



X- + 



pres 



where we used the fact that —y/2f3/{A'K) ^ a < and we posed a = —\a\. 

We use the following notation for the real and imaginary part of the resonance of H^, El 

be — i'je- From theorem [s] we have that for e small enough —l3 <he < (3 and \he — {l3 — (47ra;)^) | < 

Ce^, moreover < 7^ < Ce^, for some positive constant C . 

As we are interested on the decay-rate of a x+-state of the spin we will examine the evolution of 

the x+ component of the state at any time. Let / be the interval centered in ii^o,+ of amplitude 

2A, e^ <^ 2A ^ (47ra;)^ . We shall analyze the evolution of the following projector 
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^7++ t;x,a;' := j^^ / ^dX . 18 

647r4 X X' A |(^A^^/(47r2) - |a|) (v/A+^/(47r2) - |a|) - £2p 



Since AgJ, /5 — A>0 the integrand function in (|18j) is bounded. In turn this imphes that the 
operator ,^/++(t) is Hilbert-Schmidt for all t > 0. 



Using the substitution ^ = ^ |a| we can rewrite (18) as 



(4^)2|a|(|x|+|x'|) r J(^)g-*(/3-(47r)2(5+|a|)2)tg_(47r)25{|x| + |a;'|) 

[t; x,x ) = —jr III, I / . ^. , ^1 TT^ 2fl ^o ;rTi 17 ^'^^ 



-e^e- 



where 

/(O = v/2/3-(47r)2(e+|a|)2 (^ + |a|) 



and J' 



v/(47r|a|)2-A I I v^(47r|a|)2+A , , 

47r I" I' 47r I" I 



The four roots of the denominator of the in- 



tegrand function are easily analyzed. Two of them are real C,i{b) , C,2{£) ^ ^, close to the 

2/3 
(47r)2 



roots of the equation ^^ + 2|q;|^ — j^^ = and are outside the integration region. The last two 



Csi^)^Qi^) ^ ^ ^^6 complex conjugate, close to the roots of the equation j^^^'^ + 2e'^\a\^ + e'^ = 
0; the one corresponding to E^"^^ (let us say ^3(e)) has positive imaginary part and for e small 
enough 1^3(5) | ^ C e^. With the notation introduced above we can write 

_52g-(4^)2|a|(|x| + |x'|) r J(^^)g-i(/3-(4^)'«+l"l)')tg-(4^)25(N + |x'|) 

Notice that the density of states is approximately Lorentzian in / (for the relevance of a 

Lorentzian density of states close to a resonance see [25j and [26j). The difference with a pure 

Lorentzian behavior (Breit and Wigner for Physicist) is indicated by the presence of a quartic 

term in the denominator and of a slowly varying function in the numerator of the integrand 

function. 

In the next theorem we recover the standard result about the time behavior of metastable states. 

We extract the exponential term and estimate the remainder associated to the non-Lorentzian 

part of the density of states. 

In the following we shall consider the open subset of C defined byQ = {z E Cl 1^1 < 

diam(/ ) , ^ ai'g(-2) ^ tt}. The set Q is a semicircle in the upper complex plane, with 

center in the origin and having / as diameter. 



Theorem 4. Let ^i^^^{t) be defined like in (18), then 

^/,++(t):=L(t,5) + B(t,5) 
where the following notation was used 

L{x,x';t,e) = - ,, /f ^' , K{x,x')e-^''^'e-^^' 

with 



e-(^")'H(N+k'l) _ \a\^/2(3 - {47i\a\y 

i\ yX, X ) I j-j -j , p ——. , , r 

\x\\x'\ ^i{e)^2{£) 

and the following estimate holds true 

\\B{t,e)\\Hs<Ce\ 
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Proof. The integrand function appearing in (19) can be analytically continued to Q, and by 
the residues theorem we get 



V+(i):=^i(i)-^2(i) 



where 



Pi(t; x,x') 



^^2g-ib,tg-7,t 



n{Ue)-Cs{e)) 



K{x,x' 






and 



P2{t;x,x') 



271- 



;K{x,x') 



f{z) e'' (/3-(4'r)"(2+|a|)2) t ^-{i7T)^z{\x\ + \x'\) 



-dz . 



Where we used ,^3 = a//5 — b^ + ije — \a\. Obviously the operators -Pi(t) and P2it) depend on 

e even if it is not indicated explicitly. 

We analyze first the term P2{t,x,x'). For all x, x' G M.^ and for e small enough 

\P2it,x,x')\ ^ £2^|ir(x,a;')|e-('")''°f(^')(l"l+l"'l^ 

Where we used the fact that for all z G dQ\I , Im ((-2 + |a|)^) > 0, and that for e small enough 

1/(^)1 ^ C and \{z - U^)){z - Ue)){z - Ue)){z - QieW ^ C . 

Since |a| + inf(/ ) > we get the estimate ||-P2(^)||_h"5 ^ Ce^ . 
From the definition of L{t,e) and -Pi(t) we have 



\L{t, e; X, x') — Pi{t; x, x') 



^^2g-i6,tg-7.t 



K{x,x') 



|a|v/2/5-(47r|a|)2 fiU^)) e-(4")'«3(e)(N+l-'l) 



The estimate \\L{t,e) — Pi{t)\\Hs ^ Ce'^ comes directly from \^3{e)\ ^ Ce"^. 

The proof of the theorem is then obtained by setting B{t,e) = Pi{t) — L{t,e) — P2{t)- 



D 



As pointed out by many authors (see, e.g., [15] and references therein) the estimate given in 
theorem |4] does not make explicit the short time behavior of the solution. Notice in fact that 
the error term B{t,e) might be in principle much larger than \L{0,e) — L{t,e)\ for < t < 1. 
The following representation of ^/,++(t) is more suitable to examine details of the short time 
evolution of the projector. 



Proposition 1. Let ^/^_|_+(t) be defined like in (18) and ^ t ^ 1/e, then 

^i,++it; X, x') = a{x, x') - b{x, x'){l - e"* ^^'" *) + c{x, x')eH + 0{eH^)d{x, x') 



(20) 



where 



b{x, x') = 



27r2 



:K(x,x') 



fiOe 



^(4^)2e(|x|+|x'i) 



d^ 



/(^3(5))e-(4'^)'«3(^)(l^l+l^'l) 



27^HUe) - Qie)) 



K{x, x' 



iUe) - UeMiUe) - Ue)) 



c(x, x') 



27r2 



K{x,x') 



f{z) I {f3 - [Anfiz + |a|)2) e-(4-)'-(l-l+l-'l) 



Moreover 



WWhs ^ C ; 



dQ\l' 

HS ^ C ] 



z^ + 2\a\z^ - j^z^ - 2e^\a\z - e^ 



dz . 



WWhs^C; 



HS 



^C. 



(21) 
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Proof. The proof is obtained by writing ,^7^++(t) as ,^/,++(0) + (^/,++(t) — ^/_++(0)). Then 
a = ^/,++(0) and 

Following what was done in the proof of theorem |4] we can use the analytic continuation theory 



and the residues theorem. Then the term b in (20) corresponds to the contribution from the 
resonant pole whereas the third term is the derivative in t = of the integral on the semicircular 
part of the boundary of Q. 



The estimates (21) can be easily obtained by an argument similar to the one used in theorem 



(111), for this reason we omit the details. D 

Notice that none of the integral kernels in the representation of <^/.++(t) given in the proposition 
correspond to a projection operator. 
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